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,^ ■ Abstract 

I ' We study the complete 0{a(.w) electroweak (EW) corrections to the production 

^~^. of three Z'^-bosons in the the framework of the standard model(SM) at the ILC. The 

(-H ' leading order and the EW next-to- leading order corrected cross sections are presented, 

and their dependence on the colliding energy y/s and Higgs-boson mass nxa is ana- 
lyzed. We investigate also the LO and one-loop EW corrected distributions of the 
transverse momentum of final Z" boson, and the invariant mass of Z^Z^-pair. Our 
Q"^ ' numerical results show that the EW one-loop correction generally suppresses the tree- 

\^ , level cross section, and the relative correction with mn = 120 GeV{lhO GeV) varies 

^ ■ between -15.8%(-13.9%) and -7.5%(-6.2%) when ^i goes up from 350 GeV to 

'^. \ 1 TeV. 
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I. Introduction 

To discover the signature of new physics beyond the standard model(SM)[Tl [2] is one of 
the main goals for the forthcoming colhder experiments. The precision measurements 
of the trihnear gauge-boson couphngs are helpful for verification of non-abelian gauge 
structure, and the investigation of the quartic gauge-boson couplings can either confirm 
the symmetry breaking mechanism or present the direct test on the new physics beyond 
the SM[3J. The direct study of quartic gauge-boson couplings requires the investigations 
of the processes involving at least three external gauge-bosons. In Refs.[41 El El El [S], the 
precise predictions for the VVV productions at hadron colliders were provided. It shows 
that the QCD corrections increase the WW Z cross section at the LHC by more than 70%, 
and the QCD corrections to ZZZ production at the LHC increase the LO cross section 
by about 50% [6j. Thus, any quantitative measurement of the concerned gauge couplings 
will have to take QCD corrections into account. 

Due to heavy backgrounds, the precise measurement at a hadron collider is more 
difficult than at linear collider. The proposed International Linear Collider (ILC) by the 
particle physics community will be built with the entire colliding energy in the range of 
200 GeV < ^/s < 500 GeV and an integrated luminosity of around 500 (fb)^^ in four 
years. The machine should be upgradeable to ^/s ~ 1 TeV with an integrated luminosity 
of 1 (ab)^^ in three years[9j. Among all the ILC physics goals, the verification of gauge 
theory in the SM and finding the evidence of new physics via experimental measurement 
of the electroweak gauge boson couplings are crucial tasks too. The measurement will be 
able to be improved considerably at ILC compared with at Fermilab Tevatron and CERN 
LHC, and therefore a precise understand of the SM phenomenology at ILC to at least 
one-loop order is necessary[Tn]. Without the accurate theoretical predictions and reliable 
error estimates for important observables at ILC, it is impossible to interpret experimental 
data properly. 

The process of ZZZ production with the subsequential leptonic decays of vector bosons 



at ILC is not only an important process as a background for various new physics processes, 
but also possible to provide further tests for the quadrilinear gauge boson couplings, 
including the four gauge boson coupling, such as between ZZZZ^ which does not exist at 
tree-level in the SM, because this kind of quadrilinear couplings would induce deviations 
from the SM predicted observables|llj. 

In this paper we present the calculations of the cross sections for the process e^e~ — > 
Z^Z^Z^ at the leading order(LO) and involving complete electroweak (EW) one-loop 
(©(aeu,)) corrections. The paper is organized as follows: In the next section we present 
the calculation descriptions for the tree-level process e^e'' -^ Z^Z^Z^ . The calculation 
of the electroweak corrections at one-loop level is reported in section III. The numerical 
results and discussions are given in section IV. In the last section we give a short summary. 

II. Leading-order e+e -^ Z^Z^Z^ process 

In the calculations of the tree- level and one-loop level cross sections for the e+e~ — > 
Z^Z^Z^ process, we use the 't Hooft-Feynman gauge. The analytically calculation of the 
leading order cross section for e+e^ -^ Z^Z^Z^ process is presented in this section. We 
describe the lowest order e^e~ -^ Z^Z^Z^ process adopted for evaluating the cross section 
as 

e+(pi) + e-fe) ^ ^°(P3) + ^°(P4) + ^°(P5), (2.1) 

where pi (i = 1 — 5) label the four-momenta of incoming positron, electron and outgoing 
.Z^^'-bosons, respectively. Since the mass of electron/positron is negligible comparing with 
the colliding energy and the Yukawa coupling strength between Higgs/Goldstone and 
fermions is proportional to the fermion mass, in our work we ignore the contributions of 
the Feynman diagrams involving the couplings of H^ — e'^ — e~ and G^ — e~^ — e~ . We 
depict the tree- level Feynman diagrams contributing to the cross section of the production 
process of three Z'^-bosons at the ILC in Fig(TJ There we have 9 generic tree-level diagrams 
for the process e~^e~ -^ Z^Z^Z^ in the framework of the SM. All the Born-level diagrams 




Figure 1: The generic tree-level Feynman diagrams for e^ e -^ Z Z Z process 



can be grouped in two different topologies. Figsma-c) belong to s-channel, Figs{T][d-i) are 
grouped in t(u)-channel. The differential cross section for the process e"'"e~ -^ Z^Z^Z^ at 
the tree-level is then obtained as 



d^LO = ^\Y.\^ 



LO\ 



'd^. 



(2.2) 
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where M.lo is the amplitude of all the tree-level diagrams in Fig(TJ The factors ^ and 
J are due to three identical final Z^-bosons and spin-averaging of the initial particles, 
respectively. The summation in Eq. (j2.2p is taken over the spins of the initial and final 
particles, and d^^ is the three-particle phase space element defined as 
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III. Electroweak {0{aew)) corrections 

The 0{aew) order electroweak corrections to the Born-level e^e~ — > Z^Z^Z^ process 
consist of two parts, i.e., 

• The virtual contributions to the leading order process e"'"e~ — > Z^Z^Z^ from the 
electroweak one- loop and their corresponding counterterm diagrams; 



• The contribution from the real photon emission process e^e~ — > Z^Z^Z^ 7. The soft 
photon emission process e^e" -^ Z^Z^Z^ (7) consists IR singularities, which will be 
cancelled by the IR singularities in the contributions of the one-loop diagrams. There 
is no collinear IR singularity since we keep the nonzero mass of electron(positron); 

In the following subsections, we describe in detail the calculation procedure and discuss 
the calculation of each contribution part. 

III..1 Virtual corrections 

There are totally 2313 electroweak one-loop and corresponding counterterm Feynman 
diagrams being taken into account in our calculation, and they can be classified into self- 
energy, triangle, box, pentagon and counterterm diagrams. We depict some representative 
samples among 66 pentagon diagrams in Figl2j In our calculation of the electroweak 
{0{aew)) corrections to e^e" -^ Z^Z^Z^ process, all the one-loop Feynman diagrams and 
their relevant amplitudes are created by using FeynArts 3.3|12j. and the Feynman ampli- 
tudes are subsequently implemented by applying FormCalc5.3 programs[13j and our in- 
house routines. The electroweak one- loop amplitude involves five point tensor integrals up 
to rank 4. The numerical calculation of the integral functions(n < 4) are implemented by 
using the expressions presented in Refs.[lH[T5]. We use our independent Fortran subrou- 
tines following the expressions for the scalar and tensor five-point integrals in Ref.|16j. and 
find agreement with LoopTools2.2[13j. The Grace2.2.1 program[l7] is used to accomplish 
five-body phase-space integration for hard photon radiation process e~^e~ — > Z^Z^Z^j . 

The total unrenormalized amplitude corresponding to all the one-loop Feynman di- 
agrams contains both ultraviolet (UV) and infrared (IR) singularities. We regulate all 
singularities adopting dimensional regularization(DR) scheme |21] where the dimensions 
of spinor and space-time manifolds are extended to L> = 4 — 2e. The relevant fields are 
renormalized by taking the on-mass-shell (OMS) scheme |221 I23j. The IR singularity is 
regularized by introducing a infinitesimal fictitious mass m^. All the tensor coefficients 




Figure 2: Some representative pentagon Feynman diagrams for the process e^e 



of the one-loop integrals can be calculated by using the reduction formulae presented in 
Refs.fTOl I24j. As we expect, the UV divergence contributed by virtual one-loop diagrams 
can be cancelled by that contributed from the counterterms exactly both analytically and 
numerically. 

III. .2 Real photon emission process e'^e~ -^ Z^Z^Z^'j 

In our calculation for one-loop diagrams, there exists soft IR divergence. In order to get an 
IR finite cross section for e"'"e^ -^ Z^Z^Z^ up to the order of 0{a'^^), we should consider 
the 0{a(,y;) corrections to e^e~ -^ Z^Z^Z^ process due to real photon emission. The 
soft IR divergence in virtual photonic corrections for the process e~^e~ — > Z^Z^Z^ can be 
exactly cancelled by adding the real photonic bremsstrahlung corrections to this process 
in the soft photon limit. In the real photon emission process 

e+(pi) + e-{p2) ^ Z'ips) + Z^p,) + Z\p^) + 7(^6), (3.1) 

a real photon radiates from the electron/positron, and can be soft or hard. The general 
phase-space-slicing (PSS) method [25] is adopted to isolate the soft photon emission sin- 
gularity part in the real photon emission process e^e~ — > Z^Z^Z^^ , and the cross section 
of the real photon emission process (|3.ip is decomposed into soft and hard terms 



Aa^-eai = Act5 + Acr/f = crLo(55 +'^h)- (3.2) 

where the 'soft' and 'hard' describe the energy nature of the radiated photon. The energy 
Eq of the radiated photon in the center of mass system(c.m.s.) frame is considered soft 
if Eq < AE, and hard if Eq > AE, respectively. Then both as and an should depend 
on the arbitrary soft cutoff 6s = AE/Eh, where Ej, is the electron beam energy in the 
c.m.s. frame and equals to y/s/2, but the total cross section of the real photon emission 
process areai is cutoff AE/Eh independent. Since the soft cutoff AE/Ef, is taken to be 
a small value in our calculations, the terms of order AE/Ef, can be neglected and the 
soft contribution can be evaluated by using the soft photon approximation analytically 



[221 [231 [26] 



,A.. = -<i,„f^ / ^ (^5_ _ ^!^ I , (3,3) 



Our calculation demonstrates that the IR singularity in the soft contribution from 
Eq. (l3.3p is cancelled exactly with that from the virtual photonic corrections. Therefore, 
Acr^ + Aareah the sum of the 0{a^^^) cross section corrections from virtual, soft and 
hard photon emission contribution parts, is independent of the cutoff value 6s- The hard 
contribution, which is UV and IR finite, is computed by using the Monte Carlo technique. 
Finally, the electroweak corrected cross section for the e~^e~ -^ Z^Z^Z^ process up to the 
order of 0[a.%^) can be obtained by 

(Ttot = CTLO + Aa^, + Aareal = CfLO (1 + ^tot) ■ (3.4) 

III. .3 QED and total 0{aew) order corrections 

In analyzing the originations of the electroweak corrections, we split the full 0{aew) 
corrections to the process e^e~ — > Z^Z^Z^ into two parts, the QED correction part 
and the weak correction part. Correspondingly we define the total relative correction as 
6tot = ^QED + ^weak- The QED correction part is contributed by the diagrams with virtual 
photon in loop, and real photon emission process e^e~ — > Z^Z^Z^^ . For the countert- 
erms involved in the QED contribution, the electron/positron wave function renormaliza- 
tion constants include only photonic contribution. The remainders of the total virtual 
electroweak corrections belong to weak correction part. With above definitions we can ex- 
press the full one-loop electroweak corrected total cross section as a summation of several 
parts. 

ftot = CTLO + ^<^v,QED + Ads + AcT/i + Aa^^weak 

= CTLO (1 + 6qeD + Suieak) = (^LO (1 + ^tot) , (3.5) 

where Acr^ and Aag are the cross section corrections contributed by the virtual electroweak 
one-loop diagrams and the soft photon emission process respectively, Aa^^QED, Ao"s, Aa^ 



and Aa^^yj^ak ^-re the corrections from the virtual QED contribution, the soft photon 
emission process, the hard photon emission process and the virtual weak contribution, 
separately. Sqed, ^weak and 5tot are the relative corrections contributed by the QED 
correction part, the weak correction part and the total electroweak correction, respectively. 

As we mentioned above, there exist both ultraviolet (UV) divergency and infrared(IR) 
soft singularity in the contributions of the electroweak one-loop diagrams for e^e~ — > 
Z^Z^Z^ process, but no collinear IR singularity since we keep the nonzero mass of elec- 
tron/positron in our calculation of one-loop order corrections. After doing the renor- 
malization procedure, we verified that the UV singularity is vanished, and the IR soft 
divergency appeared in the virtual correction is cancelled by the the soft photon emission 
process e^e~ — > Z^Z^Z^-^ . 

In order to discuss the origin of the large correction when the colliding energy is 
close to the threshold of the production of three Z^-bosons, we discuss the photonic 
(QED) corrections and the genuine total electroweak corrections separately. The QED 
corrections comprise two parts: the QED virtual corrections Aa^^QED which contributed 
by the loop diagrams with virtual photon exchange in loop and the corresponding QED 
parts of the counterterms, and the real photon emission corrections Aareal- Therefore, the 
QED relative correction Sqed can be expressed as 

SqED = Sv,QED + ^real, (3.6) 

where 8v,qed = ^(yv^QEol^yLO-, and the genuine weak relative correction b^ can be got 
from 



Ow = Otot — OqeD- 

IV. Numerical results and discussion 

For the numerical calculation, we take the input parameters as follows [27] : 

me = 0.51099892 MeV, m. = 105.658369 MeV, rrir = 1776.99 MeV, 



(3.7) 



ruu = 66 MeV, m^ = 1.25 GeV, mt = 174.2 GeV, 

ma = 66 MeV, m^ = 95 MeV, m;, = 4.7 GeV, 

mw = 80.403 GeV, mz = 91.1876 GeV. (4.1) 

There we use the experimental value of W-boson mass as input parameter, but not the 
mw evaluated from G^ as in Uew scheme |28j. We take the electric charge defined in the 
Thomson limit aew(O) = 1/137.036 and the effective values of the light quark masses (m^ 
and rriii) which can reproduce the hadronic contribution to the shift in the fine structure 
constant aew{fn\) [5S]. As we know that the LEP II experiments provide the lower 
limit on the SM Higgs mass as 114.4 GeV at the 95% confidence level from the results 
of direct searches for e+e~ — > Z^H^ production [5Ul |3I], and the electroweak precision 
measurements indicate indirectly the upper bound as mn ^ 182 GeV at the 95% C.L., 
when the lower limit on niH is used in determination of this upper limit |31j. Therefore, 
in our numerical evaluation it is reasonable to take the mass of Higgs-boson being in the 
range of 115 GeV < rriH < 170 GeV. Then we shall not encounter the resonance problem 
of Higgs-boson during our calculation. 

We checked the correctness of the numerical results of the LO cross section for process 
e+e" -^ Z°Z°Z° , byusingGrace2.2.1[17] and FeynArts3.3/FormCalc5.3 ^^ packages 
separately. In adopting Grace2.2.1 and FeynArts3.3/FormCalc5.3 programs, we used both 
't Hooft-Feynman and unitary gauges separately in the calculations of the LO cross section 
to check the gauge invariance, and we got coincident numerical results. The numerical 
results of the LO cross section for the process e~^e~ — > Z^Z^Z^ , by using 't Hooft-Feynman 
gauge and taking y^ = 500 GeV and the other input parameters shown in Eqs. (|4.ip . are 
listed in Table [TJ There it is shown that there is a good agreement between the numerical 
results by adopting different packages. 

In the one- loop calculation, we must set the values of the IR regulator m^, the fictitious 
photon mass, and soft cutoff 5s = AE/E), besides the parameters mentioned in Eqs. (|4.ip . 
As we know, the total cross section should have no relation with these two parameters if 

10 



mniGeV) 


o-Lo(/^) (Grace) 


(^Loifb) (FeynArts) 


115 


1.0056(4) 


1.0055(2) 


120 


1.0139(4) 


1.0138(2) 


150 


1.0975(4) 


1.0975(2) 


170 


1.2565(4) 


1.2564(2) 



Table 1: The numerical results of the LO cross sections for the process e^e~ — > Z^Z^Z^ 
by using Grace2.2.1 and FeynArts3.3/FormCalc5.3 packages separately, and taking the 
input parameters as shown in Eqs. (j4.ip and ^/s = 500 GeV. 



the IR divergency does really vanish. Our numerical results show that the cross section 
correction at 0{a^^) order Aatot = Acrreal + '-^(^v is invariable within the calculation errors 
when rriH = 120 GeV, y/s = 500 GeV , 5s = 10^^ and the fictitious photon mass m^ varies 
from 10-1^ GeV to lO^^ GeV. 

FiglS^a) presents a verification of the correctness of our calculation for process e"'"e~ -^ 
Z^Z^Z^ including 0{aew) order corrections. The amplified curve for /S.atot including 
calculation errors is depicted in Figl3][|b). Both figures are to show the independence of the 
total ©(oew) electroweak correction on the soft cutoff 5s, when we take mn = 120 GeV 
and ^/s = 500 GeV. From Figl3l[|b) we can say that the total EW relative correction 
Aatot has no relation to the value of 5s within the calculation error range. In the further 
calculations, we set m-y = 10^^ GeV and 5s = 10^^. 

In Table [21 we list some representative numerical results of the LO and one-loop 
EW corrected cross sections (ctlo, atot), the QED and total EW corrections to the cross 
sections for e+e" — > Z^Z^Z^ process(A(7Q£;i3, Aatot), and their corresponding relative 
corrections((5Q£;D, 5tot) when ^/s = 500 GeV and the values of Higgs-boson mass are 
taken to be 115 GeV, 150 GeV and 170 GeV separately. From these data we can see that 
the one-loop EW corrections suppress the LO cross section of the process e^e~ — > Z^Z^Z^ 
, and the relative corrections are about minus few percent. We can conclude also the LO 
and EW corrected cross sections increase with the increment of mn, while the absolute 
total EW correction decreases when the value of Higgs-boson mass goes up. 
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Figure 3: (a) The dependence of the 0{acw) correction to cross section of e^e~ -^ Z^Z^Z^ 
on the soft cutoff 5s with mn = 120 GeV and \/s = 500 GeV. (b) the amphfied plot of 
the curve for the total correction Aatot in Fig-3(a), where it includes calculation errors. 



mH{GeV) 


^LoUh) 


(^totifb) 


Ao-QEoifb) 


Aatotifb) 


Sqed{%) 


Stot{%) 


115 


1.0055(2) 


0.9159(7) 


-0.0451(7) 


-0.0896(7) 


-4.49(7) 


-8.91(7) 


150 


1.0975(2) 


1.0194(8) 


-0.0444(8) 


-0.0780(8) 


-4.04(7) 


-7.11(7) 


170 


1.2564(2) 


1.1989(9) 


-0.0393(8) 


-0.0575(9) 


-3.12(7) 


-4.58(7) 



Table 2: The numerical results of cr^o, atot, ^ctqed, Ao"tot(in femto bar), and their 
corresponding relative EW and QED corrections(5tot5 ^qed) for the process e+e~ -^ 
Z^Z^Z^ , when ^ = 500 GeV and mu = 115 GeV, 150 GeV, 170 GeV respectively. 
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Figure 4: (a) The LO(o"io)) 0{aew) EW, QED corrected cross sections((Ttoi) ctqed) for 
the process e~^e~ — > Z^Z^Z^ as the functions of cohiding energy ^/s with mn = 120 GeV, 
150 GeV separately, (b) The corresponding relative EW, QED relative corrections ((^tofj 
Sqed) versus ^/s. 

The numerical results of the LO, 0{aew) EW, QED corrected cross sections (o-^^Oi crtot, 
o'qed) and the total relative EW, QED corrections (^tot, Sqed) for the process e~^e~ — > 
Z^Z^Z^ with rriH = 120 GeV, 150 GeV as the functions of colliding energy ^/s are 
plotted in FigsHJa) and (b) respectively. As indicated in Figl^a), The curves for the cross 
sections of (Tloi (^tot and ctqed increase quickly in the ^/s region of [350 GeV, 550 GeV] 
and decrease when ^/s > 600 GeV . The two figures show the 0{aew) corrections always 
suppress the corresponding LO cross sections of process e~^e~ -^ Z^Z^Z^ in both cases 
of niH = 120 GeV and mn = 150 GeV separately, but the QED correction parts can 
enhance the LO cross sections when ^/s > 700 GeV. We can read out from Figlljb) that 
the corresponding total EW relative corrections for mn = 120 GeV and 150 GeV vary in 
the ranges of [—15.8%, — 7.5%] and [—13.9%, — 6.2%] respectively, when ^/s runs from 
350 GeV to 1 TeV. We can see also from these two plots that in the colliding energy ^/s 
region near the threshold of the production of three Z'^-bosons, the main contribution to 
the total EW relative correction(5tot) comes from the QED correction part. That is due to 
the Coulomb singularity effect coming from the instantaneous photon exchange in loops 
which has a small spatial momentum. 
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Figure 5: The distributions of the transverse momenta of Z'^-bosons (P^) at the LO and 
up to EW one-loop order with mn = 120 GeV and y/s = 500 GeV . 

We present the distributions of the transverse momenta of final Z^-bosons at leading 
order and up to one-loop order, daio/dp^, da^Lo/dp^', when mn = 120 GeV and ^/s = 
500 GeV in FigjSl There we pick the p^ of each of the three Z'^-bosons as an entry in this 
histograms, then the final result of the differential cross section is obtained by multiplying 
factor 1/3. In this figure we can see that the EW one-loop correction suppresses obviously 
the LO differential cross section daio/dp'^ when p^ > 50 GeV , but the EW correction is 
small when p^ < 25 GeV . It also shows that the EW corrections do not observably change 
the LO distribution line-shape oip^, and both the differential cross sections of daio/dp^ 
and daNLo/dp^ have their maximal values at about p^ ~ 50 GeV respectively. 

We plot the distributions of the invariant mass of Z^Z^-pair, denoted as Mzz, at the 
LO and up to EW one-loop order with tuh = 120 GeV and ^/s = 500 GeV in FigEl 
Here we can see that the EW correction slightly enhances the LO differential cross section 
when Mzz < 250 GeV, but suppresses daLo/dMzz obviously when Mzz > 250 GeV. 
The suppression of daio/dMzz is due to the contribution from the hard photon emission 
process e^e~ -^ Z^Z^Z^j at the 0{a^^) order, in which the momentum balance between 



14 









1 ' 


1 ■ 1 


1 


1 


. 




0.005- 




JP 


ijisiirC'^^^^f^ 


;;^^ 




- 


5 


0.004 - 




/^ 




''■■■'■">-. 


\ 




1. 






/ 






:... 


, \ 


■ 






fj 


0.003 - 




r 


LO 






M 


- 


S 


NLO 


T3 


0.002 - 




: 


m„=120GeV 
*=500GeV 






■) 


- 




0.001 - 












- 




n nnn 












" 




u.uuu - 

i; 


>o 


1 ' 1 ' 1 

200 250 300 


1 

350 


1 

400 


45 



M„ [GeV] 



Figure 6: The distributions of the invariant mass of Z^ Z^-paii^Mzz) at the LO and up 
to EW one-loop order when mn = 120 GeV and a/s = 500 GeV. 

the sum of the momenta of three Z°-bosons and that of the radiated hard photon will 
reduce the value of invariant mass Mzz and change the line-shape in the range with large 
Mzz. 

V. Summary 

In this paper we describe the impact of the complete one-loop EW corrections to the 
scattering process e^e~ -^ Z^Z^Z^ in the SM. This channel can be used to measure the 
quartic vector boson coupling at ILC. We investigate the dependence of the LO, 0{aew) 
EW and QED corrected cross sections on colliding energy ^/s and Higgs-boson mass, and 
present the LO and EW one-loop corrected distributions of the transverse momenta of 
final Z'^-bosons and the LO and EW corrected differential cross sections of invariant mass 
of Z^Z^-pair. To see the origin of some of the large corrections clearly, we calculate the 
QED and genuine weak corrections separately. We conclude that both the Born cross 
section and the EW corrected cross section for e+e~ -^ Z^Z^Z^ process are sensitive to 
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the Higgs boson mass in the range of 115 GeV < mn < 170 GeV. We find the 0{aew) 
corrections generally suppress the LO cross section, the LO distribution of the momenta of 
Z^-bosons and the LO differential cross sections of invariant mass of Z^Z^-paii for process 
e^e~ — > Z^Z^Z^ . Our numerical results show that when mfj = 120 GeV{150 GeV) and 
the colliding energy goes up from 350 GeV to 1 TeV, the relative EW correction varies 
from -15.8%(-13.9%) to -7.5%(-6.2%). 
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